Spiral ground state against ferroelectricity in the frustrated magnet BiMnFe206 
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The spiral magnetic structure and underlying spin lattice of BiMnFe2 06 are investigated by low- 
temperature neutron powder diffraction and density functional theory band structure calculations. 
In spite of the random distribution of the Mn^^ and Fe^^ cations, this centrosymmetric compound 
undergoes a transition into an incommensurate antiferromagnetically ordered state below Tat — 
220 K. The magnetic structure is characterized by the propagation vector k = [0, P, 0] with j3 ~ 
0.14 and the P22i2il'(0/30)0s0s magnetic superspace symmetry. It comprises antiferromagnetic 
helixes propagating along the 6-axis. The magnetic moments lie in the ac plane and rotate about 
7r(l + /3) ~ 204.8 deg angle between the adjacent magnetic atoms along b. The spiral magnetic 
structure arises from the peculiar frustrated arrangement of exchange couplings in the ab plane. The 
antiferromagnetic coupling along the c-axis cancels the possible electric polarization, and prevents 
ferroelectricity in BiMnFe2 06. 

PACS numbers: 75.25-j, 61.05.fm, 61.66. Fn, 75.30.Et 



I. INTRODUCTION 

The coupling between magnetism and ferroelectric- 
ity is one of the intriguing phenomena in solid-state 
physics.'^ Apart from the ongoing studies of the underly- 
ing microscopic mechanisms ,1213 the effect itself is relevant 
for applications,^ and stimulates extensive experimental 
work on diverse systems varying from bulk transition- 
metal compound^^m ^ heterostructures.'^ Magnetoelec- 
tric effects in bulk systems typically conform to one of 
the two following scenarios: (i) the magnetism arises from 
transition-metal cations with a partially filled d shell, and 
the ferroelectricity is driven by lone-pair cations, such as 
Bi'^+ or Se''+ (Refs. [5] and [H]), or (ii) electronic effects [a 
spiral (helicoidal) magnetic structure or a charge order- 
ing] break the symmetry and cause ferroelectricity.'^Ii^'^ 
The former mechanism ensures large electric polariza- 
tion, which, however, is weakly coupled to the mag- 
netism. The latter scenario provides a strong coupling, 
but a small electric polarization. The combination of the 
two approaches is clearly advantageous, but difficult to 
achieve. The best known example is BiFeOa,-^ which 
exhibits a plethora of interesting effects relate d to t he 
coupling between ferroelectricity and magnetism.'^^Ei! 

The recently discovered complex oxide BiMnFe2 06 
represents another system combining the two scenarios 
plausible for multiferroicity: the lone-pair Bi'^'^ cation 
and the spiral magnetic ground state. Neither of the 
two, however, lead to ferroelectricity. The polar dis- 
placements of Bi^+ are ordered in an antiferroelectric 
mannerjli^ whereas the spiral magnetic structure is non- 
polar due to a strong antiferromagnetic coupling along 
the crystallographic c direction. 




FIG. 1. (Color onhne) The crystal structure of BiMnFe206. 
The Fe and Mn cations are situated in the oxygen octahedra 
(in the color version, blue and violet polyhedra denote the 
FeMnl and FeMn2 positions, respectively). 

The unique crystal structure of BiMnFe2 06 features 
fragments of the hypothetical hep oxygen-based MO 
building blocks that are related by a mirror opera- 
tion into a polysynthetically twinned structure.!^ The 
Fe'^"'" (d^) and Mn^+ (c?^) cations are nearly randomly 
distributed over two crystallographically distinct posi- 
tions, both octahedrally coordinated by oxygen atoms 
(Fig. [T]). The octahedra are interconnected into a frame- 
work through corner-, edge-, and face-sharing, provid- 
ing complex paths for the magnetic exchange. In spite 
of the random distribution of the Fe and Mn cations, a 
long-range magnetic order was reported in BiMnFe206 
below T/v — 220 K, according to ^^Fe Mossbauer spec- 
troscopy, magnetic susceptibility, and heat capacity mea- 
surements. The low-temperature neutron powder diffrac- 
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FIG. 2. (Color online) Neutron diffraction patterns for 
BiMnFe206 at different temperatures. The arrows show the 
magnetic reflections. 



tion revealed the incommensurate propagation vector of 
the magnetic structure, yet no details of the ground-state 
spin arrangement have been reported.!^ 

Spiral magnetic structures are capable of inducing 
ferroelectricity in a number of transition-metal ox- 
ides, such as TbMn03_^ NiaVaOs,'^^ MnW04,^ and 
Bao.sSro.sZnaFeiaOaa.'^ In contrast to the aforemen- 
tioned compounds, BiMnFeaOg remains paraelectric be- 
low Tjv.l^ To understand the lack of ferroelectricity, we 
studied the magnetic structure and explored the underly- 
ing frustrated spin lattice. Although the strong frustra- 
tion leads to the formation of magnetic helices propagat- 
ing along b, the unfrustrated antiferromagnetic (AFM) 
coupling along c induces the 222 point symmetry of the 
magnetically ordered state and cancels the possible elec- 
tric polarization. 



II. METHODS 

The powder sample of BiMnFeaOe was synthesized by 
high-temperature solid-state reaction in air. Stoichiomet- 
ric amounts of raw materials (BiaOa, MnaOa, FcaOa) 
were ground thoroughly and heated up to 800 °C in 
5 h. After annealing for 10 h, the powder was reground, 
pressed into a pellet, and heated at 1000 °C for 100 h 
with several intermediate grindings. 

Neutron powder diffraction (NPD) data were collected 
on a — 12 g sample of BiMnFcaOg, contained in a 9.5 mm 
diameter vanadium can. A closed cycle He refrigerator 
was used for temperature control. Patterns were col- 
lected with the BT-1 32-detector high-resolution neu- 
tron powder diffractometer at the National Institute of 
Standards and Technology Center for Neutron Research, 



Gaithersburg, MD. A Cu(311) monochromator with a 
90° takeoff angle and 15 min in-pile coUimation was used. 
The neutron wavelength was 1.5402(1) A. Data from the 
32 detectors were combined to give pseudo one-detector 
data over a total scan range of 3° < 26* < 167.75° with 
a step size of 0.05° {29). The magnetic structure was 
analyzed with the JANA2006 program.'^!' 

Density functional theory (DFT) band structure calcu- 
lations were performed in the full-potential local-orbital 
(FPLD) code.'^Si We used local density approximation 
(LDA)P supplied with a mean-field (DFT+U) correction 
for correlation effects in the Fe/Mn 3d shell. The fc-mesh 
comprised 192 points for the 40-atom unit cell and 64 
points for the 80-atom supercell. The correlated shell 
was parameterized with an effective on-site Cou lomb re- 
pulsion Ud = 5 eV and exchange Jd = 1 eVj ^^ ^ ^^ l whereas 
the double-counting was corrected in the fully-localized- 
limit (atomic limit) fashion. To evaluate individual ex- 
change couplings, total energies for a number of collinear 
spin configurations were mapped onto a classical Heisen- 
berg model. The validity of the computational results 
was checked by calculations within generalized gradient 
approximation (GGA^PH and by choosing Ud values of 
4 and 6 eV. Similar to Cu^"*" oxides,'^^! the exchange- 
correlation potential (LDA vs. GGA) has marginal effect 
on the spin model, whereas a change in Ud only shifts 
the exchange couplings in a systematic way and keeps 
their ratios nearly constant. The DFT-based spin model 
was further studied by classical Monte-Carlo simulations 
with the spinmc algorithm of the ALPS package.'^S' 



III. RESULTS 
A. Magnetic structure 

At room temperature (RT), BiMnFeaOe crystallizes 
in an orthorhombic unit cell with a = 5.03590(3) A, 
6 = 7.07342(4) A, c = 12.65425(6) A, and Pbcm space 
symmetry.f^^l Below T/v — 220 K, extra reflections ap- 
pear on the NPD patterns (Fig. 2]). These reflections are 
of magnetic origin and can not be attributed to a struc- 
tural phase transition because the x-ray diffraction exper- 
iment does not show any change down to T = 120 K.^^^ 
The magnetic reflections on the T — 7 K pattern are 
indexed with a propagation vector k = [0, /3,0] with 
/3 — 0.1379(1). The propagation vector is inside of the 
Brillouin zone and has a star with two arms k, — k. The 
little group of the propagation vector Gk is Pb2im. In 
BiMnFcaOg, there are two symmetrically independent 
magnetic species in the nuclear Pbcm structure: FeMnl 
(8e: 0.4891,-0.1597,-0.6036) and FeMn2 (4a: 0,0, i) 
(Fig. [T]). In the little group Pb2im, this corresponds 
to two four-fold magnetic sites for the FeMnl position 
and one four-fold magnetic site for the FeMn2 position 
(origin at 0,0,^). There are four one-dimensional irre- 
ducible representations (irreps) for the propagation vec- 
tor k = [0, /3, 0] in the space group Pbcm; their char- 
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TABLE I. Characters of the irreducible representations of the little group of the propagation vector 
k = [0,/3,0] for the space group Pbcm {a — e^^^) and the corresponding magnetic superspace 
groups. The decomposition of the magnetic representation Fmag for all four-fold magnetic sites is: 
bmag = Sbi + 3r2 + SFg + 3r4. 
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TABLE II. Symmetry operators for the P&cml'(0/30)sOOs and Pfecml'(0/30)00ss magnetic superspace groups for origin shift 
5 — Q,\ along 14. The operators resulting from combination with the (1' | 000 |) generator, always present, are not shown 
(— m means a "time inversion" operation, while m is an operation without time inversion). The operators forming the common 
subgroups are printed in boldface. 
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acters are given in Table |Tj The magnetic structure is 
transformed according to one of the irreps or their com- 
bination. Alternatively, we can describe the symmetry 
of the incommensurately modulated magnetic structure 
by embedding it into a higher-dimensional space and 
applying magnetic superspace groups defined in (3-1-1)- 
dimensional superspace.'^SHSSl 'pj-^g magnetic moment on 
atom i is expressed as a vector function: 

N 

Mi{x4) = Mio+^[Mi„sSin(27rna:4)+Mi„cCOs(27rna;4)], 

n=l 

where n denotes terms of the Fourier series, 
X4 = k(T -I- Ti) is an internal coordinate, T is the 
lattice translation of the nuclear structure, and 
is the position of the atom i in the unit cell of the 
nuclear structure. Monitoring of the intensity of the 
magnetic and nuclear reflections upon varying temper- 
ature revealed that there is no magnetic impact into 
the intensity of the nuclear reflections, and therefore 
Mio = for all magnetic sites. Since only the first-order 
satellites were observed, the Fourier series are reduced 
to the n = 1 terms. A magnetic (Shubnikov) superspace 
group describing the transformations of the magnetic 
modulation waves can be set in correspondence with 



each irrep. The magnetic superspace groups are based 
on generators of the little group Gk, but the symmetry 
elements transforming the propagation vector k to — k 
also should be taken into account. This yields four 
possible magnetic superspace groups listed in Table [l] 
The explanation of the magnetic superspace group 
symbols is provided in Ref. [33] 

All four magnetic superspace groups were tested 
in the refinement. Acceptable solutions were found 
in magnetic superspace groups PbcmV {Of30)sOOs and 
Pbcml' {0(30)00ss, both with the same reliability factor 
for magnetic reflections R^^^ — 0.057. In spite of the 
relatively low reliability factor, the correspondence be- 
tween the experimental and calculated NPD profile was 
far from ideal for both models. This indicates that the 
actual solution requires a combination of the irreps 
and and can be realized in a common subgroup of 
P6cto1'(0/30)sOOs and P6cto1'(0/30)00ss. This subgroup 
depends on a relative shift S along the internal space of 
the conventional origins of two superspace groups. The 
list of operators of P6cml'(0/30)sOOs and their intersec- 
tions with the operators of P&cml'(0/30)00ss for S = 0,\ 
are provided in Table |ll] The resulting common sub- 
groups are P12i/cll'(0/30)sOs and P22i2il'(0/30)0s0s. 
For any 5 values inequivalent to the cases mentioned 
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above the common subgroup is P12ill'(0/30)ss. 

The solutions in the P12i/cll'(0/30)sOs and 
P22i2il'(0/30)0s0s groups provide the same quahty 
of the Rietveld fit and can not be distinguished on this 
basis. The P12i/cll'(0/30)sOs magnetic superspace 
symmetry results in a collinear magnetic structure with 
an antiferromagnetic transverse amplitude modulated 
wave with the magnetic moments confined to the 
ac plane. The amplitude of the magnetic moment 
modulation varies from almost zero to 5.6 /is, which 
is unrealistically high for the Mn'^+ and Fe'^"'" cations. 
Thus, the P12i/cll'(0/30)sOs solution was ruled out. In 
the P22i2il'(0/30)0s0s model, there are three magnetic 
symmetrically unequivalcnt atoms: FeMnla, FeMnlb, 
and FeMn2 (see Table III I, all at the general four-fold 
sites of the P22i2i space group. The (l'|000 i) 
operator of the magnetic superspace group requires 



the Mio term to be zero and constrains the magnetic 
moment modulation functions to odd harmonics only, 
thereby resulting in the absence of the magnetic impact 
into intensity of the nuclear reflections and the absence 
of even-order magnetic satellites.'^'' The components 
of the magnetic moment modulation function for the 
general four-fold site are related by symmetry elements 
of the P22i2il'(0/30)0s0s magnetic superspace group, 
as shown in Table IIVI 

Since in the nuclear structure the atoms FeMnla and 
FeMnlb are crystallographically equivalent, it is reason- 
able to assume that the modulation of the magnetic mo- 
ment at these positions follows the same type of mod- 
ulation waves. The solution was found with the follow- 
ing restrictions (not imposed by the magnetic superspace 
group) on the coefficients of the magnetic moment mod- 
ulation functions: 



M^,^ (FeMnla) = Mc,:^ (FeMnla) = 

M^,:r (FeMnla) = -Mc,^ (FeMnla) = 

Ms,^(FeMn2) Mc,^(FeMn2) 

Mc,2(FeMn2) = M,,^(FeMn2) 



M^,^ (FeMnlb) = M^,^ (FeMnlb) = Mi 
M^,:r (FeMnlb) = Afc,2 (FeMnlb) = M2 





The refined magnetic moment components along the b 
axis were smaller than their standard deviations and were 
fixed to zero. The refined magnetic parameters at dif- 
ferent temperatures are provided in Table [V] The ex- 
perimental, calculated, and difference NPD profiles at 
T = 7 K are shown in Fig. [3] The temperature depen- 
dencies of the ordered magnetic moment for the FeMnl 
and FeMn2 positions are given in Fig. |4j Extrapolations 
of these dependencies with the M = Afo[(l - T/TnTY 
function give coinciding Tn values of 221(2) K and 
219(4) K for the FeMnl and FeMn2 positions, respec- 
tively, and the magnetic moments of A//o (FeMnl) = 
3.76(4) and A/o(FcMn2) = 4.0(1) /^s. 

The spin arrangement in the BiMnFe206 magnetic 
structure is shown in Fig. [5] BiMnFc206 adopts a spi- 
ral magnetic structure consisting of antiferromagnetic 
helixes propagating along the 6-axis with a period of 
~ 3.56. The helixes are associated with the FeMnl 
and FeMn2 atomic chains running along the 6-axis, 
where the magnetic atoms are separated by h/2. The 
magnetic moments rotate for 7r(l -I- /3) ~ 204.8 deg 
about the 6- axis on going between the adjacent mag- 
netic atoms in the chains. The refined parameters A^i, 
Af2, and represent the magnetic moment of the 



FeMnl atom 



, the magnetic 



A4^ (FeMnl) 

moment of the FeMn2 atom [Af(FeMn2) = jAfa]], and 
the phase shift between the FeMnl and FeMn2 helixes 
(tan(p = -^). 
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FIG. 3. (Color online) Experimental, calculated, and differ- 
ence NPD pattern for BiMnFe2 06 at T = 7 K. The black 
(dark) and green (light) bars mark the reflection positions for 
the nuclear and magnetic structures, respectively. 



The refined structural parameters and interatomic dis- 
tances at r = 7 K (space group P22 i2i ) an d RT (space 
group Pbcm) are listed in Tables III and VI respectively. 
The crystal structure at T = 7 K was refined with fixed 
parameters of the magnetic structure. The crystal struc- 
tures at both temperatures are virtually identical indi- 
cating that the magnetic ordering does not infiuence the 
nuclear structure. With the P22i2il'(0/30)0s0s magnetic 
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TABLE III. Crystallographic parameters and atomic coordinates in BiMnFe206 at T = 7 K and room temperature (RT). 
Atomic displacement parameters (7iso are given in 10^^ A^. The symmetry operators of the P22i2i space group are listed in 
Table HI] 







T = 


7 K 






RT 




Space group 




P22 


i2i 






Pbcm 




a (A) 




5.02305(7) 






5.03589(3) 




MA) 




7.06232(9) 






7.07341(4) 




c(A) 




12.6424(2) 






12.65425(6) 




Atom 


X 


y 


z 


(7ieo 


X 


y z 




Bi 


0.9711(3) 


-0.1301(2) 


0.751(1) 


0.27(4) 


0.9705(1) 


-0.1305(1) 1 


0.68(2) 


FeMnlcj^] 


0.486(3) 


-0.160(2) 


-0.607(1) 


0.32(8) 


0.4891(3) 


-0.1597(2) -0.60355(7) 


0.80(3) 


FeMnlb'" 


-0.488(3) 


0.160(3) 


0.607(1) 


0.32(8) 








FeMn2'' 


0.000(6) 


0.000(6) 


0.500(2) 


0.32(8) 





\ 


0.37(4) 
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0.157(3) 


0.439(3) - 


-0.6349(9) 


0.47(2) 
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0.4372(3) -0.6366(2) 
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0.86(7) 
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0.3444(6) 
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FeMnla and FeMnlb: 0.7Fe+0.3Mn; FeMn2: 0.6Fe+0.4Mn 



TABLE IV. Symmetry-imposed relations between the components of the magnetic moment modu- 
lation functions for the general four-fold site of the P22i2i l'(0/30)0s0s magnetic superspace group 
with the origin at 0,j,0,0. 
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FIG. 4. (Color online) Temperature dependence of the mag- 
netic moment for the FeMnl and FeMn2 positions. 



symmetry, the space group of the average nuclear struc- 
ture should be P22i2i (point group 222). Although the 
spiral magnetic ordering eliminates the inversion center, 
it does not create a polar direction. Indeed, no indication 
of ferroelectricity below Tn was found in BiMnFc2 06 by 
dielectric permittivity measurements.^ 



B. Electronic structure 

Owing to the complex crystal structure, an em- 
pirical assignment of individual exchange couplings in 
BiMnFe206 is a formidable challenge. The problem can 
be solved by electronic structure calculations that eval- 
uate individual couplings and, therefore, establish a reli- 
able microscopic magnetic model. Recent studies prove 
the remarkable accuracy of DFT for diverse correlated 
systems, including frustrated magnets with highly intri- 
cate spin lattices.l^SESl Prior to computing exchange inte- 
grals, we discuss the electronic structure of BiMnFe206, 
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TABLE V. Refined magnetic parameters for BiMnFe206 at difi'erent temperatures (see text for 
notations). Reliability factors are listed in the order of Ri (overall), Ri (nuclear reflections), Ri 
(magnetic satellites), Rp. 
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7 


3.25(2) 


1.94(2) 


3.98(3) 


3.80(3) 


0.13801(8) 


0.017, 0.014, 0.028, 0.036 


57 


3.13(2) 
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3.71(2) 


3.65(3) 


0.13740(4) 


0.016, 0.013, 0.025, 0.030 


100 


3.09(4) 


1.75(3) 


3.50(4) 


3.55(5) 


0.1316(1) 


0.025, 0.022, 0.040, 0.048 


150 


2.74(6) 


1.40(5) 


2.57(6) 
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0.1232(2) 
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200 


1.69(6) 


0.88(5) 


1.56(6) 


1.91(8) 


0.1175(3) 


0.029, 0.027, 0.044, 0.051 


215 


0.95(9) 


0.54(9) 


0.68(9) 


1.1(1) 


0.1155(9) 


0.028, 0.028, 0.036, 0.049 



TABLE VL Selected interatomic distances (in A) in BiMnFe206 at T = 7 K and room temperature (RT). 
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1.99(2) 










FeMnlb-02b 


1.99(2), 2.51(2) 
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FeMnlb-04 


1.92(2) 






RT 












Bi-01(x2) 


2.207(2) 


FeMnl-Ol 


1.928(3) 


FeMn2-01(x2) 


1.964(2) 


Bi-01(x2) 


2.686(2) 


FeMnl-02 


1.973(3) 


FeMn2-02(x2) 


2.056(2) 


Bi-04 


2.246(3) 


FeMnl-02 


1.992(3) 


FeMn2-03(x2) 


2.057(2) 


Bi-02(x2) 


2.696(2) 


FeMnl-04 


2.007(2) 






Bi-04 


2.816(3) 


FeMnl-03 


2.029(3) 










FeMnl-02 


2.489(2) 







which is a key to understanding the magnetic behavior. 

Computational analysis of BiMnFe206 is complicated 
by the intrinsic disorder of Fe and Mn atoms. Unfortu- 
nately, state-of-the-art computational tools, such as vir- 
tual crystal approximation (VGA) or coherent potential 
approximation (CPA), cannot be applied to this case, 
since they provide an averaged description of a disor- 
dered system, while we are targeting magnetic proper- 
ties that depend on local interactions. The insulating 
nature of BiMnFe2 0j^ suggests localized 3d electrons 
of transition metals; therefore, the distinct spin-| Fe'^+ 
(d^) and spin-2 Mn'^+ (d^) sites are randomly distributed 
on the spin lattice. Further, orbital ordering for Mn'^"'" 
results in dramatic differences among the Fe-0~Fe, Mn- 
0-Mn, and Fe-O-Mn superexchanges (see below). An 
averaged VCA or CPA picture would not reproduce any 



FeMn2 




't tt^ < •* 

V V V V t 



FIG. 5. (Color online) The arrangement of magnetic mo- 
ments in BiMnFe2 06. The unit cell for the nuclear structure 
is outlined. 
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FIG. 6. (Color online) LDA density of states for BiFesOe- 
The Fermi level is at zero energy. The gapless energy spec- 
trum is caused by underestimation of the electronic correla- 
tions in LDA. 



of these features. To capture effects arising from the 
localized moments of Fe^+ and Mn^+ , we calculated ex- 
change integrals for several systems with ordered Fe and 
Mn atoms, and accessed all possible scenarios of the 
superexchange. In these calculations, we kept the ex- 
perimental atomic positions, but imposed different ar- 
rangements of Fe and Mn, spanning purely Fe (BiFesOg) 
and purely Mn (BiMnaOg) cases as well as intermediate 
(BiMnFe206 and BiMn2Fe06) configurations. 

The LDA band structure depends only marginally on 
the Fe/Mn ordering. The density of states (DOS) fea- 
tures Bi 6s bands below —10 eV, O 2p states between 
—7 eV and —2 eV, transition- metal 2>d states at the Fermi 
level, and Bi 6p states above 3 eV (Fig. [6]) . The effect of 
substituting Mn for Fe is a change in the electron count 
and the ensuing shift of the Fermi level within the id 
bands. Irrespective of the Fe/Mn ratio, LDA band struc- 
tures are metallic due to the heavy underestimation of 
strong electronic correlations. The insulating spectrum 
is correctly reproduced by DFT-|-f7 (see upper panel of 
Fig. [t]). In particular, the band gaps of about 1.5 eV 
for BiFegOe and 1.1 eV for BiMugOg at [7^ = 5 eV are 
consistent with the black color of the compound. 

Superexchange couplings in insulating compounds in- 
timately depend on the orbital state of the transition 
metal. Five unpaired electrons of Fe'^+ fill five d orbitals 
and leave no orbital degrees of freedom. By contrast, 
Mn'^"'" has four unpaired electrons only, and hence one 
of the d orbitals is unoccupied in the Mott-insulating 
state. Since there is a sizable crystal-field splitting (about 
1.5 eV) driven by the octahedral coordination of the 
FeMnl and FeMn2 sites, the three t2g states are half- 
occupied, while the two Cg states have to share the re- 
maining electron. The half-filled orbital is picked out by 
a weak distortion of the octahedral local environment. In 
the following, we analyze such distortions in more detail 
to determine which of the two eg orbitals is half-filled in 




-4 
Energy (eV) 

FIG. 7. (Color online) Top: LSDA-H(7 DOS for one spin 
channel in the ground-state AFM configuration of BiMnaOg 
{Ud = 5 eV). Middle and bottom: orbital-resolved DOS for 
3d states of Mn^^. Four out of five d orbitals show filled 
spin-up and empty spin-down states, whereas the states of 
the remaining d orbital {d^2_y2 for Mnl and d^^2_^2 for Mn2) 
are mostly above the Fermi level for both spin directions. The 
Fermi level is at zero energy. 



the insulating state. 

The FeMnl position re vea ls one long bond of about 
2.5 A (FeMnl-02, Table |vi]) that we choose as the lo- 
cal z axis. According to simple electrostatic arguments 



of crystal-field theory, the long bond along z shifts the 
d^;,2_j,2 orbital down in energy with respect to the (i^2_j,2 
orbital. The opposite scenario is found for FeMn2, where 
the local environment is a squeezed octahe dron . The 
short FeMn2-01 bond (about 1.96 A, Table |vi|, which 



we take as the local z-axis, shifts the d3j,2_^2 orbital up in 
energy. Therefore, the unpaired electron of Mn'^"'" should 
occupy the d^^^^r^ orbital for FeMnl and the dx2-y2 or- 
bital for FeMn2. This conclusion is verified by DFT-I-C/ 
calculations that place unpaired electrons on the respec- 
tive orbitals, while shifting both spin-up and spin-down 
states for the remaining Cg orbital {d^^-y^ and d^z^^r^ for 
FeMnl and FeMn2, respectively) above the Fermi level 
(see Fig.[7f.37 

The orbital order in BiMngOe is formally of the ferro- 
type within the FeMnl and FeMn2 sublattices, yet of 
the antiferro-type between the two sublattices. However, 
this notation is deceptive because of the different local 
axes on neighboring atoms. The magnetic model largely 
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deviates from the conventional scenaricl^ of AFM su- 
perexchange for ferro-type orbital order and ferromag- 
netic (FM) superexchange for antiferro-type orbital or- 
der. In the following, we show that the peculiar arrange- 
ment of empty eg orbitals in BiMnaOg induces sizable 
FM superexchange for nearly all couplings and alters the 
spiral ground state that arises from the purely AFM spin 
lattice of BiFesOe- 



C. Microscopic magnetic model 

Exchange couplings calculated for different superex- 
change scenarios (Fe-O-Fe, Mn-O-Mn, Fe-O-Mn, and 
Mn-O-Fe) are listed in Table |Vl!| The Fe-O-Fe, Mn- 
0-Mn, and Fe-O-Mn superexchanges show sharp dif- 
ferences for most of the couplings, while the Fe-O-Mn 



and Mn-O-Fe cases are only different for interactions 
between the FeMnl and FeMn2 sublattices (i.e., when 
the two metal sites are not related by symmetry). In 
Table VII we restrict ourselves to short-range couplings 
matching direct connections between the FeMn octahe- 
dra. Long-range interactions are expected to be weak, as 
confirmed by the following qualitative argument. Long- 
range superexchange requires suitable overlap of atomic 
orbitals along a M-O-O-M (or even more complex) path- 
way, and can be achieved for one orbital channel only. By 
contrast, short-range superexchange is possible for most 
of the orbital channels, and should therefore dominate 
in systems with several magnetic orbitals. To verify this 
conclusion for BiMnFe2 06, we calculated the exchange 
couplings within the crystallographic unit cell and within 
a supercell doubled along a. The resulting values of Ji 
agreed within 10% and indicated weak long-range cou- 
plings in BiMnFe206. Below, we demonstrate that our 
minimum microscopic model, restricted to short-range 
couplings, is sufficient to explain the spiral magnetic 
structure of BiMnFe206 and the lack of ferroelectricity 
in this compound. 

The spin lattice of BiMnFe2 06 incorporates nine in- 
equivalent exchanges (Fig. [s] and Table VII I. Eight of 
these couplings are found in the ah plane, whereas Jg 
connects the layers along c. The interlayer coupling is 
weakly influenced by the Fe/Mn substitution, and re- 
mains one of the leading AFM interactions for all su- 
perexchange scenarios (Table VII). The robustness of J% 
should be traced back to the orbital order for Mn'^+ . The 
local z-axis of the FeMnl octahedra roughly matches the 
crystallographic c direction. Therefore, the replacement 
of Fe by Mn results in an empty dx2-y2 orbital (Fig. [t]) 
that contributes weakly to the superexchange along c. 
By contrast, most of the intralayer couplings are heavily 
affected by changing Mn for Fe. 

In BiFeaOe, eight intralayer couphngs are AFM. Ex- 
cluding the apparently weak Ji , we arrive at seven AFM 




O FeMnl 
• FeMn2 



FIG. 8. (Color online) Spin lattice (left) and the respec- 
tive part of the crystal structure (right) showing frustrated 
exchange couplings in the ab plane. Open and filled circles 
denote the FeMnl and FeMn2 positions, respectively. 



released by Mn'^+ that renders most of the intralayer cou- 
plings FM, and reduces the remaining AFM couplings 
below 13 K (Table [yn]). The Fe-O-Mn and Mn-O-Fe 



couplings ranging from 19 K to 74 K (Table VII). Tri- 
angular loops abound (Fig. [8| and lead to the strong 
frustration of the spin lattice. This frustration is largely 



cases are intermediate, with partially reduced AFM cou- 
plings. Large FM contributions to the superexchange can 
be traced back to empty d orbitals of Mn'^^. These or- 
bitals provide the strong cr-overlap with oxygen orbitals 
and a leading contribution to the superexchange. 

To evaluate the ground state of the proposed model, 
we performed classical Monte-Carlo simulations for a 
16 X 16 X 16 finite lattice with periodic boundary con- 
ditions. The unit cell of the spin lattice incorporated six 
magnetic atoms within one layer (half of the crystallo- 
graphic unit cell). The temperature was set to 10 K, well 
below the ordering temperature (see further in this 
section). Spin-spin correlation functions for purely Fe'^^ 
{S = I) an d purely Mn^+ {S = 2) lattices are given in 
Table IVIIll The correlation functions normalized for 
are —1 for antiparallel spins and -1-1 for parallel spins. In- 
termediate values indicate non-collinear configurations. 

We consider the Fe'^"'" case first. The normalized spin- 
spin correlation for Jg is close to —1; therefore, the inter- 
layer ordering is collinear and AFM. Similar correlations 
on the J5 and Jq bonds indicate the same propagation 
vector along b for the FeMnl and FeMn2 sublattices. Fur- 
ther on, similar correlations on the J2 and J3 as well as 
on J7 and Jg bonds signify the same magnetic order along 
the respective bonds and, consequently, a twice shorter 
periodicity along b for the magnetic unit cell compared 
to the crystallographic unit cell. To find out the ordering 
pattern along a, we note that the respective FeMn2 atoms 
are connected via J2 and J3 or J7 and Jg bonds. The cor- 
relations along these paths are different; therefore, the 
spins on the FeMn2 atoms should be parallel (i.e., the 
magnetic moment rotates for a certain angle ipi on the 
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TABLE VII. Interatomic distances (in A) and leading exchange couplings (in K) in BiMnFe206 
calculated with the supercell procedure (LSDA+(7, Ud ~ 5 eV) for diflerent scenarios of superex- 
change. Negative Ji denotes FM coupling. The intralayer couplings Ji — J7 and Jg are depicted in 
Fig. [8] whereas Jg connects the layers along c. 





Distance 




Exchange couplinf 


;s Ji 






Fe-O-Fe 


Mn-O-Mn 


Fe-O-Mn/Mn-O-Fe 


Jl 


2.916 


1 


-1 


-50 


J2 


3.010 


19 


-9 


29/10 


J3 


3.101 


25 


-38 


2/-23 


Ja 


3.462 


26 


7 


-38 


Js 


3.537 (FeMn2) 


30 


13 


43 


Je 


3.538 (FeMnl) 


61 


-8 


48 


Jl 


3.685 


37 


-1 


23/30 


Jg 


3.706 


74 


79 


127 


Jg 


3.759 


66 


-12 


20/70 



strong enough to induce the spiral order. While BiFeaOg 
TABLE VIII. Normalized spin-spin correlations (5i5,)/5 BiMnaOg present two opposite scenarios of strong 

BiFea Oa (Fe^^, g ^ f ) and BiMnaOa (Mn^+, g = 2). ^^^^ frustration, respectively, BiMnFe^Og hes be 

tween these distinct regimes. The spiral ground state 
of this compound is driven by the frustration of in- 
tralayer exchange couphngs that remain predominantly 
AFM for m ixed Fe-O-Mn superexchange pathways (see 
Table [VIl|. 

For an additional test of our microscopic model, we 
calculated magnetic susceptibility (x) and estimated the 
Neel temperature T/v a-s the position of the kink in the 
temperature dependence of x- Transition temperatures 
for BiFegOe (420 K) and BiMngOg (150 K) reasonably 
agree with the experimental Tjv of 220 K that lies be- 
tween the two calculated values. 





BiFeaOe 


BiMnsOe 


Jl 


0.288 


0.957 


J2 


0.240 


0.954 


J3 


0.234 


0.974 


J4 


-0.835 


0.945 


Js 


-0.693 


0.940 


Je 


-0.706 


0.951 


Jl 


-0.827 


0.945 


Js 


-0.992 


-0.993 


Jg 


-0.835 


0.956 



J2 bond and for the opposite angle —(pi on the J3 bond). 
Thus, the propagation vector of the magnetic structure 
is [0, /3, 0], in agreement with the experiment. The prin- 
cipal spin arrangement is described by three parameters, 
which are the angles between the spins on the J5 (Jg), 
and J J bonds. According to our simulations, these angles 
are (f = 226 deg, ipi — 76 deg, and (p2 = 146 deg, respec- 
tively, in remarkable agreement with the experimental 
values oi if — 206 deg, ipi = 39 deg, and = 167 deg. 
Our microscopic model for BiFcgOg reproduces the ex- 
perimental magnetic structure of BiMnFe206 quite well. 
The remaining discrepancies are likely related to the par- 
tial replacement of Fe by Mn. 

The complete substitution of Mn for Fe changes the 
magnetic ground state. According to Table VIII most 



of the normalized spin-spin correlations in BiMngOg are 
close to -1-1. The only negative correlation refers to the 
Jg bond and indicates AFM interlayer coupling. The 
intralayer ordering is now collinear FM due to predom- 
inantly FM exchange couplings (Ji, J2, J3, Jg, J7, and 
Jg, see Table |VII[ ). The spin lattice is still frustrated 
by AFM couplings J4 and J5, which, however, are not 



IV. DISCUSSION AND SUMMARY 

Our experimental and computational study provides 
microscopic insight into the physics of BiMnFe206. This 
compound features a strongly frustrated spin lattice with 
predominantly AFM exchange couplings that induce the 
spiral magnetic order. The comparison between BiFeaOg 
and BiMngOg suggests that AFM exchange couplings are 
an essential prerequisite for the spiral magnetic order. 
The crystal structure itself may still allow for different 
ground states; therefore, it is the transition-metal cation 
that determines the type of the long-range magnetic or- 
der. In this respect, first reports on cation substitutiorP^ 
look promising, because the change in the Fe/Mn ratio 
or an incorporation of other magnetic cations could alter 
the magnetic structure and other physical properties. 

BiMnFe206 is a magnetic compound that conforms 
to two mechanisms of ferroelectricity. First, lone pairs 
of Bi^"*" induce polar displacements. Second, the non- 
coUinear magnetic structure breaks the inversion symme- 
try and allows for magnetic-field-induced electric polar- 
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ization. However, none of the two mechanisms succeed in 
rendering the compound ferroelectric. Bi'^^-related po- 
lar displacements form an antiferroelectric pattern that 
leads to zero net polarization. The electronic mechanism 
meets a similar obstacle of the AFM interlayer coupling 
and the ensuing non-polar (albeit non-centrosymmetric) 
magnetic structure. The lack of polarity in both atomic 
and magnetic structures naturally explains the absence 
of the ferroelectric response in BiMnFe2 06 below T^]^ 

Our results demonstrate that the simple criteria of fer- 
roelectricity and multiferroicity are not universal, since 
the polarization induced by any kind of polar distor- 
tion (atomic displacement or spin spiral) can be wiped 
out by an overall antiferroelectric/antiferromagnetic or- 
der. The combination of lone-pair and transition-metal 
cations does not necessarily lead to a magnetic ferro- 
electric, whereas an arbitrary incommensurate magnetic 
structure may not allow for the electronic mechanism 
of ferroelectricity. These simple observations make the 
search for multiferroics a challenge, and put forward the 
charge-ordering mechanism as a more robust approach to 
the design of magnetoelectric materials. 

The cancellation of polarity in an incommensurate 
magnetic structure has been proposed for the spin-chain 
cuprate NaCu20j^ and fo r a n umber of layered com- 
pounds, such as a- CaCraOj^QP^and a-SrCr204,'^ that 
were considered as potential multiferroics. The non- 
frustrated AFM interlayer coupling seems to be a gen- 
eral obstacle for ferroelectricity in layered systems, like 
a-CaCr204,'*'^ or in three-dimensional systems with two- 
dimensional frustrated units, as in BiMnFe206. To over- 
come this problem, one has to design materials with 



FM or frustrated interlayer couplings. In the case of 
BiMnFe206, cation substitutions preserve the long-range 
ordei^^ and can be promising for tuning this compound 
toward ferroelectric and possibly multiferroic behavior. 

In summary, we have solved the spiral magnetic struc- 
ture of BiMnFe2 06 and proposed a microscopic magnetic 
model for this compound. The ground state features 
spins lying in the ac plane and propagating along 6 by a 
~ 204.8 deg rotation about the b axis. The two inequiva- 
lent FeMn positions reveal the same propagation vector. 
The spiral magnetic structure is driven by the strong 
frustration of antiferromagnetic exchange couplings on 
a complex spin lattice in the ab plane. However, the 
coupling along the c direction is non-frustrated, leading 
to antiparallel spin arrangement in neighboring layers. 
The resulting magnetic structure is non-polar and, alike 
the centrosymmetric (antiferroelectric) atomic structure, 
precludes the ferroelectric behavior of BiMnFe206. 
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